Interacting boson problems are QMA-hard 
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Computing the ground-state energy of interacting electron (fermion) problems has recently been 
shown to be hard for QMA, a quantum analogue of the complexity class NP. Fermionic problems are 
usually hard, a phenomenon widely attributed to the so-called sign problem occurring in Quantum 
Monte Carlo simulations. The corresponding bosonic problems are, according to conventional wis- 
dom, tractable. Here, we discuss the complexity of interacting boson problems and show that they 
are also QMA-hard. In addition, we show that the bosonic version of the so-called JV-representability 
problem is QMA-complete, as hard as its fermionic version. As a consequence, these problems are 
unlikely to have efficient quantum algorithms. 

PACS numbers: 03.67.a, 05.30.Jp, 89.70.Eg 



Many important model Hamiltonians in physics, such 
as the Hubbard model (both fermionic and bosonic) and 
those for superconductivity and the quantum Hall effect, 
involve at most two-body interactions The ground- 
state wavefunction and energy of these Hamiltonians play 
a key role in understanding these fascinating phenom- 
ena. In some of these phenomena, electrons are the ma- 
jor players. Problems involving fermionic particles such 
as electrons often seem to be computationally more dif- 
ficult than those with their bosonic counterparts. This 
intractability is often attributed to the so-called "fermion 
sign problem" , occurring in Quantum Monte Carlo sim- 
ulations . On the other hand, bosonic problems do not 
suffer from the sign problem [3[ and are thus regarded as 
tractable. 

Indeed, Schuch and Verstraete and Liu, Christandl 
and Verstraete have recently shown that computing the 
ground-state energy of general interacting electrons is 
QMA-hard The complexity class QMA (Quantum 

Merlin- Arthur) is a generalization of the class NP (non- 
deterministic polynomial time) to the quantum realm. It 
was introduced by Kitaev @ in the study of the so-called 
local Hamiltonian problem, where, roughly speaking, the 
goal is to determine the ground-state energy of a spin 
Hamiltonian involving only few-body interaction terms. 
QMA-hard problems are regarded as difficult, unlikely to 
be solved efficiently even by a quantum computer. How- 
ever, a quantum computer, if given the solution to any 
problem in QMA, along with a suitable "certificate" or 
"witness state", can efficiently verify whether the solu- 
tion is correct or not. In fact, as a result of a series 
of works 0, H, Hi , even for nearest-neighbor two-body in- 
teractions among spin- 1/2 particles on a two-dimensional 
lattice, the local Hamiltonian problem is QMA-complete. 
With higher magnitude of spins in one dimension, the 
local Hamiltonian problem can be QMA-complete as 
well [To| . Understanding and classifying the complexity 



of physical models and investigating hard problems using 
statistical mechanical tools have become important re- 
search endeavors 0, 1, 1, 0, !, @M El El EI 0, EE EH , 
as a result of interplay between physics, mathematics and 
computer science. 

The fact that interacting fermion problems are hard 
motivates us to investigate the corresponding bosonic 
problems. The complexity of bosonic problems seem 
less explored. Could bosonic problems be so hard as to 
be intractable? We show that generic nearest-neighbor 
two-body interacting boson problems are indeed QMA- 
hard. Inspired by the recent work of Liu, Christandl, 
and Verstraete [5( that shows QMA-completeness of the 
(fermionic) A^-representability problem [T3, EH , we study 
its bosonic version. In an A-represcntability problem one 
is given a two-particle reduced density matrix p and needs 
to decide whether there exists a global A^-body state a 
that is consistent with p. The JV-representability prob- 
lem has been vastly explored in quantum chemistry [l9| , 
and its solution would enable efficient solution of ground- 
state energy for generic two-body interacting fermionic 
systems. We show that the bosonic A-representability 
problem is also QMA-complete. Similarly, we show that 
the bosonic A-rcprcscntability problem given only diag- 
onal elements is NP-hard. 

QMA-hardness of interacting boson problems. Let us 
consider boson creation and annihilation operators , a,j 
for the j'th site or mode, and denote by |f2) the vacuum 
state without any bosons. The creation and annihilation 
operators obey the following commutation relations [2(| : 

[a,, a,-] = = [<4,aj], [di, at] = (1) 

The use of these operators preserves the symmetry of the 
bosonic wavefunctions under permutations, and any N- 
boson wavefunction can be represented in the following 
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second-quantization formalism: 

iv>> = E ^...^(4r.,.(aif-\n), (2) 

U + ...+i m =JV 

where ifc (0 < ik < AT) denotes the number of bosons at 
the k th site, and m is the total number of sites. Note 
that we restrict ourselves to the subspace of states with 
exactly N bosons. 

We construct a bosonic Hamiltonian 7ibosc, whose 
ground-state energy determines the ground-state energy 
of the following quantum spin glass model TL: 

where i, j label the site, = 1 denotes the identity 
and en 1 ) = o x , cr' 2 ' = a y , and cr' 3 ) = er z are the three 
Pauli matrices. The coefficients are real but arbi- 
trary. Oliveira and Terhal Q showed that determining 
the ground-state energy of TL is QMA-hard, even if the in- 
teractions are restricted to nearest neighbor sites (i, j) in 
the two dimensional square lattice. By way of reduction, 
solving the ground-state energy of TLbosc is also QMA- 
hard. To construct 7iboso , we use the Schwinger boson 
correspondence between qubit and boson states (see, e.g., 
Ref. [21|) given by the following map: 

erf «-> ajb, + b\a t , of «-> i(&ja, - a-fcj), erf <-> aja^ - 6-^, 

. (4) 

where the operators a,i,bi correspond to distinct sites. 
(This is similar to the mapping in the fermionic case [f|.) 
It is easy to verify that the bosonic operators obey the 
commutation relations of the corresponding Pauli opera- 
tors. We can regard the qubit at site i as a single boson 
that can be in one of two different degrees of freedom: 
\zi) <-> (alY-ibjy-^in) with z t £ {0,1}. This corre- 
sponds to the dual-rail encoding of a photonic qubit in 
the Knill-Laflamme-Milburn linear-optics quantum com- 
putation scheme Hence, TV qubits can be repre- 
sented by N bosons in 2N sites (or N sites with each bo- 
son possessing two distinct internal degrees of freedom) : 

. . • , z N ) «- (a\r (bl) 1 - 21 ■ ■ ■ (ajv) 2N (^jv) 1 " 2 " I") ■ 

.(5) 

Next, we consider the resulting bosonic Hamiltonian. 
As any two-local qubit Hamiltonian can be written as a 
linear combination of terms with at most two Pauli op- 
erators, the corresponding bosonic Hamiltonian can be 
written as a combination of products of at most two anni- 
hilation and two creation operators. To ensure that there 
be exactly one boson on the pair of sites corresponding 
to i, we add the following extra terms: 

Pi = (aloi + btbi -I) 2 , (6) 

which commute with other terms in the Hamiltonian. 
The total bosonic Hamiltonian is then 

H hose = H(a^b^a,b) + J2 cP h (7) 

i 



which involves at most nearest-neighbor two-body inter- 
actions. By making the weight c of these projectors large 
enough, e.g., Yi,j,^u \ c ij\ A r (A 7 -l)/2, we guarantee that 
the ground state of the full Hamiltonian has exactly one 
boson per site. Thus TLbose may be represented with 
at most a polynomially larger number of bits as com- 
pared to TL. Thus, if one can compute the ground-state 
energy of general bosonic Hamiltonians with at most 
two-body interactions, one can solve general spin-1/2 
two-local Hamiltonian problems. As solving the latter 
is QMA-hard, solving the former is QMA-hard as well. 
This shows that interacting boson problems are generally 
difficult. 

QMA-hardness of bosonic N -representability problem. 
We consider the number of bosons N to be fixed, and 
assume that the number of modes m that the bosons oc- 
cupy is large enough, i.e., m > 5 N for some constant 
6 > 0. The number of different ways Af m that N identi- 
cal bosons can occupy m sites is Af m = ( + ™ X )i which 
grows exponentially in N, i.e., M m ^ (8+l) N /6 N when 
N is large. Given an TV-boson state p( N > , the two-boson 
reduced density matrix is calculated by tracing out all 
but two of the bosons: p*- 2 -* = Trjv-ap^) where p( N *> is 
in general a mixture of states \ip) of the form @ with 
exactly N bosons. A precise definition of p^ is given 
via its matrix elements: 

Pijll = N (N-l) ^ ajai) ' (8) 

where the bracket indicates the expectation value over 
the state p^ . Note that the one-boson reduced density 
matrix p^> = Trjy-ip^, defined via p\£ = (a\a,i)/N, is 
completely determined once p^ is known: 

I 

Informally, the bosonic iV-representability problem 
(with m sites) asks whether there is an iV-boson state 
whose two-particle reduced density matrix equals a given 
state p. To deal with technical issues related to precision, 
we are promised that when there is no TV-boson state con- 
sistent with it, every two-particle reduced density matrix 
is "far away" from p. Formally, we are given a two-boson 
density matrix p of size [m(m + l)/2] x [m(m + l)/2], and 
a real number > l/poly(7V), with all numbers specified 
with poly (TV) bits of precision. We would like to decide 
if: ("YES" case) There exists an TV-boson state cr such 
that Tr W -2 (cr) = p, or if ("NO" case) For all iV-boson 
states cr, ||TrAr_ 2 (cr) — p\\i > f3. 

We show that the bosonic N- representability is QMA- 
hard under Turing reductions [24j . In other words, we 
show that given an efficient algorithm for bosonic N- 
representability, we can efficiently determine the ground 
state energy of Hbosc, a QMA-hard problem as estab- 
lished above. In the sequel, we refer to an algorithm for 
iV-representability as the "membership oracle" . 
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The first step is to write the two-particle interacting 
terms in 7ibose m terms of a complete orthonormal set, 
Q, of two-particle observables: iitwo-body = Z^QeQ 7Q Q, 
where the number of elements I = \Q\ ~ 0(m 4 ). Note 
that poly(m) is poly(A), and so is poly(Z). The observ- 
ables Q are constructed as in the fermionic case @. We 
define aj = a i2 a il , for all pairs / = (11,12), * 1 < i<i ■ (Note 
that in the case of fermions i\ < 12-) We fix a total order 
(denoted by -<) on pairs of indices /. The observables in 
<2 are defined as follows: 



1 fat a j + a\ai), forl^J, (10) 



(OjCLj — CLjCLl), for I -< J, (11) 

(12) 



— a\ai, for / -< L 



ri[ 



where the factor rc.j = 1 if i\ ^ 12, fii = 2 if %\ = and L 
denotes the last pair in the ordering. These operators are 
Hermitian and the Xjj and Yjj have expectation values 
in the interval [— 1, 1] and the Zi have expectation values 
in [0, 1] for any two-particle state. In the two-particle 
Hilbert space, they are orthogonal to each other under 
the trace operator, e.g., Tt(XijZk) = for all I,J,K. 
They also form a basis for the two-boson density matri- 
ces, i.e., any such matrix p^ may be written as 



,(2> 



= Z L + '^2 "z, (Zi - Z L ) 

+ ( a (Xi.,) X U + » ( 13 ) 

I-<J 

where a Q = Tr(Qp( 2 )) for all Q E Q. Using Eq. © 
we see that the expectation value {a\a^) of the one- 
body terms, can be expressed as linear combination of 
the a Q . Thus we have Tr{H hoae p^) = Sqgq 7q a Q, 
where 7q includes the contribution from both one-body 
and two-body terms, and aQ are the coefficients of 
p( 2 ) = Ttn-2(p^) ■ Finding the ground-state energy 
is equivalent to minimizing the linear function f(d) = 
^q^qIqOlq subject to the constraint that a G K^, 
where Km denotes the convex set of all a such that the 
corresponding state is iV-representable. 

The above minimization of energy, subject to the con- 
vex constraint that p is iV-representable, belongs to 
a class of convex optimization problems which can be 
solved using the shallow-cut ellipsoid algorithm [25|, |26[ 
with the aid of an 7V-representability membership oracle. 
If Km is contained in a ball of radius R centred at the 
origin, and it contains a ball of radius r, the run time 
is poly( log(i?/r)) and the error in the solution due to 
computation with finite precision is l/poly(i?/r). The 
algorithm will be efficient and of polynomially bounded 
error, if R/r is at most poly(Z). 

From the discussion leading to Eq. (| 13[) it follows that 
Km is contained in a ball of radius R — \fl centered at 



the origin. However, the method used by Liu et al. [f| 
of regarding aj as creation operator for a hole in site 
j cannot be employed to suitably bound r from below. 
Instead, we explicitly construct a set of TV-boson states 
such that the convex hull of the corresponding vectors 
{oiq} contains a ball of radius l/poly(Z). Relying on 
the property that bosons can occupy the same site, for 
each Q we construct states so as to maximize or minimize 
aQ. The resulting points have the property that for any 
coordinate axis, there exist at least two points at constant 
distance along that axis. As a consequence, we show that 
their convex hull (which is contained in Km) contains a 
ball with radius r > l/poly(Z), centered at the center of 
mass of the points. 

An algorithm for bosonic iV-representability thus en- 
ables efficient calculation of the ground-state energy of 
the interacting boson Hamiltonian 7ibosc- As a conse- 
quence, iV-rcprcscntability is QMA-hard. This may come 
as a surprise, as bosonic problems are generally regarded 
as being easier than the corresponding fermionic prob- 
lems @. 

QMA-completeness. Is the bosonic iV-representability in- 
side QMA or even harder? We show that the bosonic 
A^-representability problem is indeed inside QMA, im- 
plying that the problem is QMA-complete. To establish 
this, we construct a QMA proof system, i.e., describe a 
witness state r (over polynomially many qubits) for the 
iV-representability of a given two-boson density matrix p, 
and a polynomial-time quantum algorithm V (the "veri- 
fier") that expects such a pair p,r as input. The verifier 
determines probabilistically whether a given two-boson 
density matrix p is ./V-representable, or is far from being 
iV-representable. In the "YES" case, V outputs "YES" 
with probability pi — 1, and in the "NO" case, V out- 
puts "NO" with probability > p = 1 - 1 /poly (AY The 
gap pi — po can be amplified to 1 — e -15017 ^ [q. |23|. 

For the witness state, we represent an TV-boson state 
a using m qucfits (d-dimensional quantum systems), via 
the following correspondence, a.k.a. Holstein-Primakoff 
bosons (see, e.g., Ref. (2l|): 



A,= 



1 



-St, 



1 



-St, 

where S i are the raising/lowering operators for i'th 
spin, and 2s > N. The spin operators above sat- 
isfy the bosonic commutativity relations provided the 
total spin magnitude is s. The boson number states 
|n) e {|0), |1), 1 2s)} at one site correspond to the spin 
states \s n ) <E {|s), \s — 1), | — s)}, and d = 2s + 1. 
A bosonic observable O = ajaja;afe + a\.a\ajai is trans- 
formed into 6 = A\A^AiAk+A\A\AjAi, which is a ten- 
sor product of at most four single-qucfit observables, in 
contrast to the involvement of non-local Jordan- Wigner 
string operators in the fermionic case [f|. 

The expectation value (O) can be estimated efficiently. 
One method is to explicitly diagonalize the observable O 
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as Si an d measure the given qudit representa- 

tion a of the bosonic state a in the basis Repeat- 
ing the measurement on polynomially many copies of cr, 
we can estimate (O). We note that a quantum circuit us- 
ing qudits with d = 2s + 1 can be implemented efficiently 
by an equivalent circuit using qubits. 

The witness r consists of polynomially many blocks, 
where each block has m quefits that represent a state a 
that is claimed to be an TV-boson state with TrAr_ 2 (5') = 
p. The verifier V measures, on each block, the observ- 
able J^k A^Ak = ms — J^k $k t° check whether the par- 
ticle number is N. If it is not TV, V outputs "NO". 
This measurement projects each block onto the space 
of fixed particle number states. If the particle number 
is N, V continues to perform measurements for a suit- 
able set of observables (e.g., those corresponding to Q) 
using the projected states. It compares whether the out- 
comes match the expectation values specified by p, to 
check for consistency. It outputs "YES" if the errors 
are less than /3/poly(A r ) (for a suitable polynomial), oth- 
erwise outputs "NO". When p is TV-representable, the 
prover supplies polynomially many copies of the state 
a such that Trjv_2(cr) = p, the verifier always answers 
"YES" (i.e., pi = 1). When p is not TV-representable, 
the prover can cheat by entangling different blocks of qu- 
dits. Using a Markov argument, as was first employed 



by Aharonov and Regev [27| and later by Liu (26|, one 
can show that the verifier will still output "NO" with 
probability > /3/poly(Y). Thus Y-representability is in 
QMA [H]. 

Further consequences. Now that we have shown the 
QMA-completeness of Y-representability for bosons, we 
may follow the same argument as by Liu, Christandl, 
and Verstraete Q and reach the conclusion that pure- 
state bosonic Y-representability is in QMA(fc). This is 
achieved because the essential point is to verify the purity 
of the certificate. Next, consider the problem of deciding 
bosonic TV-representability of the two-boson density ma- 
trix when only the diagonal elements = (ajajajGij) 
are specified. If one considers the case m = 2N and the 
mapping by the Schwinger representation Q, one finds 
that as in the fermionic case, the solution enables one to 
solve the ground-state energy of local spin Hamiltonians 
which only contain commuting a z operators. The lat- 
ter corresponds to a classical spin-glass problem, and is 
known to be NP-hard [l6| • Thus the problem of deciding 
iV-representability given {Dij} is also NP-hard. 
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